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Исследование стационарной конвективной 
математической модели 
кристаллизации вещества 


Исследуется стационарная задача Стефана [6 учётом конвективного движения в ЖИДКОЙ фазе на плоскости. 
Получено уравнение свободной границы в зависимости от интенсивности вихря. Строится приближенное 
решение задачи. 


Постановка задачи 


Математическое моделирование теплофизических процессов в спецметаллургии 
оперирует такими усредненными величинами, как понятие сплошной среды, её плот- 
ность, теплоемкость, температуро- и теплопроводность, вязкость и др., и, опираясь на 
законы переноса массы, импульса и энергии, составляют основные дифференциаль- 
ные уравнения для определения других усредненных характеристик — векторного поля 
скоростей, температурного поля и др. Все теплофизические параметры среды в той 
или иной степени зависят от температуры, а законы, управляющие этими зависимос- 
тями, носят экспериментально-эмпирический характер и до настоящего времени весь- 
ма неполны, особенно в области высоких температур, поэтому поневоле приходится 
рассматривать линейное приближение, применимое в случае малых перепадов темпе- 
ратуры, когда зависимостью упомянутых параметров от термодинамического состоя- 
ния можно пренебречь. Адекватность математической модели физическому явлению 
возрастает вместе с уточнением упомянутых зависимостей, вследствие чего еще большее 
значение приобретают соответствующие эксперименты общепознавательного харак- 
тера. Линеаризованная теория задач массо- и теплопереноса оперирует также с известным 
количеством безразмерных параметров изучаемого явления, таких, как числа Рейноль- 
дса, Пекле, Нуссельта, Прандтля и др., и значение этих параметров для возможно 
большего диапазона материалов и режимов тоже приобретает повышенное значение. 

Теплофизические процессы в кристаллизаторе, сопровождающиеся фазовыми пе- 
реходами вещества, описываются математической моделью, в которой температура 
каждой из фаз удовлетворяет уравнению переноса тепла со своими теплофизически- 
ми коэффициентами. На границе раздела фаз обе температуры постоянны и равны 
температуре фазового перехода (для химической однородной среды), а на заданных 
частях границы, — стенках кристаллизатора, поддоне, — поддерживается определенный 
режим (теплоотвод, теплоизоляция и др.). Поверхность раздела фаз (фронт кристал- 
лизации) является неизвестной, или «свободной» границей, и для ее определения допол- 
нительно задается «условие Стефана», означающее, что тепловой поток через фронт 
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кристаллизации в сторону твердой фазы равен тепловому потоку со стороны жидкой 
фазы плюс скрытая теплота фазового перехода. Жидкая фаза рассматриваемого про- 
цесса заслуживает специального исследования из-за априорной возможности сущест- 
вования поля скоростей, вызывающего интенсивную теплопередачу путем конвекции. 
Усиленная циркуляция в расплавленной шлаковой ванночке была обнаружена в иссле- 
дованиях академика Б.Е. Патона и его сотрудников [1]. 

Основной целью настоящей статьи является изучение гидродинамических яв- 
лений в жидкой фазе, поскольку здесь экспериментальные исследования, насколько 
известно, отсутствуют. 


Изучается стационарный случай в полосе О = {-1 <х<1 Н<у< 0}. Обозначим 
через у кривую, отделяющую жидкую фазу ру от твердой О, , при этом концы у 
лежат на вертикалях х = 1. Обе области Бу И В, предполагаются односвязными 
и симметричными относительно оси у. Пусть \(х, у) — функция тока, удовлетворяю- 
щая условиям: Аи = д, (х,у) Е В, ‚ м = сопяй > 0, и = 0,(х,у) Е дру. Здесь м счита- 
ется достаточно малым численным параметром. Требуется определить, кроме функции 


> и 
тока у(х,у), тройку (и (х,у),у) по следующим условиям: 
А. Аи“ -ууих + У хи = 0, (х, у) Е ВА = сои$й > 0; 


и* (х,0) = 0, -1<х<Ъо = сот! > 1; 


и; + @и* =0, х=+1,(х,у) е 00; ; (1) 
Аи =0, (ху) ЕО; и (х,Н)=0, -1<х<1, “ен ЕЕ, 
2 2 
Уи | —к?Ми* = 0,(х,у) Е У, 


здесь к = соп5, 0<к<1; @ —числа Нусельта. 


Приближенное решение задачи 


Предложен метод изучения нелинейной задачи (1), состоящий в разложении ре- 
шения в ряд по степеням малого параметра д [2]: 


[©.0) [©.0) 
и(х, узи) = У иКи к (х, у), и" (х.з и) = У иКин(х,у); 


к=0 к=0 


[©.0) 
Ух, и) = Хи, у : у = У(х, и), —15х <1 


к=0 
Нулевое приближение и0(х, у) ищем как минимум функционала 


Ци и ,у)= [2 +? @кдучк? [[[мз? и? | кар + 
р» р 


+к? ах [м2 1] до | [и 16, 


г в 
г, г, 


(2) 
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на соответствующем множестве допустимых функций, здесь | = дру П {х = +1 


Г, =д), П {х = +1. Функционал (2) в классе функций и >0в ру представим сле- 


дующим образом: 


2 2 Го) 
10» У,) [как [аки нод [ве - [му (Ьи)]4и+ 
А 1 А) 1 


и 2и 


+4 | (и [Би +, Си) Чи, 


где А| = (-1<х<Ъ 0<и<1), Д› = (-1<х<1 1<и<о), у1(х,и) и у2(х, и) — реше- 
ния уравнений и! (х, у) -и| =0, и>(х, у)-и_ =0. 


Будем минимизировать функционал /\ (у, У. ) при помощи сумм 


г Т, : =: р-и Ё 9, р = 
Ум акх и" +Н, (хи) ЕД}; у2и(х,и) = и. Бхи^, (х,и) Е А, 
1=0к=1 о- 1=0к=0 


применяя при этом метод Ритца [3]: 


ОГ „Б_, 
ОЧ»), Я, =0, р=1,2...,Т.;4=0,1,...,Г, 


ба 
Ва 1 г О. РОВ, 
о Е (3) 
У ЕН=о али 
к=1 к=0 ит к=1 к=0 +28 
вр аль +Я; ии | 
1=0 к= С /=0к=0 


Лемма. Пусть система Ритца (3) имеет решение при некоторых значениях параметров 
@0 =@0, 4 =@,к=к. Тогда решения этой системы ак; (@ *,@0,к), В (0,0 к) 
непрерывно зависят от параметров 4 ‚@%,к в некоторой окрестности точки (@% ‚@,к). 


Проблема сходимости приближенных решений исследована в [3], [4]. 
Теорема. Пусть // достаточно малая величина. Тогда справедлива формула: 


5) Ор), уе, (4) 
и, „(х, у) 


где у,(х) — решение уравнения и,(х, у) -1=0 в классе функций и,,(х,у)>0 вр; 


7: Уи =У@Ф-и 


и, (х, у) — решение краевой задачи: 
Ди = (х, у), (х,у)еБзи(х,0) = 0, и(х,Н) =0,0х< Е 


и, (0, у) =0,Н <у<0;и, + аи =0;х=Ъ(х, у) ед; \у; 


Ло, у) = м, У, при (х,у) ЕР,” и Лоу) =0 при (х,у)еД, . 
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Приближенный анализ влияния конвекции 
на фронт кристаллизации 

Численный анализ осуществлялся на основании формулы (4). В качестве функции 
и, (х, У) берется решение проблемы минимума функционала (2), которое может быть 


построено методом Фурье при К =1@,' =@, =@). Функции \/ (х, У) и и, (х, у) на- 
ходятся из условий минимума функционалов 
0 


(у) = [ие +2ау, ть = [ (и? и? +2 пдахау + в [ и? у), 
Бр Ур) 


Н 


на множествах 
т ={у:уеС'(Б),и =0,(х,у) 00; |, Т, ={и:ие С'(Р), и(х,0) = 0, и(х,Н)=0}. 
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Рисунок 1 — Линии кристаллизации 


Проблема минимума эффективно решается при помощи метода Ритца. При этом 
приближения Ритца \/, и и, строятся следующим образом: 


т ту 


у, = У - НХу- хп сих у, 


&=0 /= 


т Ту 


и, =У(-Н)> У ах’, п = зар (т+2т, ), 


К=0 = 0<А<т 


106 «Искусственный интеллект» 172010 


Исследование стационарной конвективной математической модели... 2 М 


где у, (х) — решение уравнения и, (х, у) -1=0, (х,у) = вклассе функций и, (х, у) > 
>0О вр. Сходимость приближенных решений к точным решениям, соответствующих 


краевых задач, изучена в [5]. Неизвестные коэффициенты а, и с, определяются из 
условий минимума функций 1 (\,) = 1 (с), Т.(и,) = Б(а,). 


Численный эксперимент осуществлялся при определенных значениях теплофи- 
зических значений параметров. На рис. 1 изображен график линий кристаллизации 


у(х, и) при различных значениях д, соответственно при д=-0,5; 0; 0,1; 0,5. Кри- 


вые у(х, и) строились в виде многочленов 


У(х, и) =а, (их +а, (и)х’ +а (и)х+а, (и). 
Вычисления производились при Н =-10, у =1,25, @, = 3,5, при этом у = у(х,0,5); 


У =у(х, 0,1); У. =у(х,-0,5) иУу-у (х) : 
Проделанный численный эксперимент подтверждает влияние конвективного теп- 
лообмена на процесс кристаллизации. Эксперимент сохранит СВОЙ смысл, если пара- 


метры @) и @, брать в некоторой малой окрестности чисел @, = 3,5, аК=1. 
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О.С. Мненко 

Дослдження стацюнарно! конвективно! математично! модел! кристалзацй речовини 
Дослджуеться стацюнарна задача Стефана з урахуванням конвективних рухв у рдиннйЙ фаз! на площин. 
Доведено равняння в1льно! гранищ в залежности в1д 1нтенсивност! вихору. Побудовано наближене рипення 
задач1. 
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Вееагсь о? Зайопагу Майетайса! Моде! о? СгузбаШтайоп 07 Зи б$апсе уу В Сопуесйоп 
Туюо-дитепз1опа| %анопагу сопуесйоп Зе п рго Мет и Налма рБазе 15 шуезизаеа. Еопла оЁ Нее Боипдагу 
едлайоп 1$ оатеч4. Те арргохипайоп зоайоп 1$ Ба|а. 
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